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Abstract 

We discuss the associated W^H T production at the CERN Large Hadron Col- 
lider. The dependence of the hadronic cross section on the Higgs sector parameters 
is investigated in detail in the framework of the general Two Higgs Doublet Model 
(THDM). We study the possible enhancement of the THDM prediction for the cross 
section compared to the prediction of the Minimal Supersymmetric Standard Model 
(MSSM). We find regions in the THDM parameter space where the THDM predic- 
tion can exceed the one of the MSSM by two orders of magnitude. These regions of 
large cross section are in agreement with theoretical bounds on the model, derived 
from the requirement of vacuum stability and perturbative unitarity, and are not 
excluded by experimental constraints. 
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1 Introduction 



Spontaneous symmetry breakdown is a necessity in theoretical descriptions of electroweak 
phenomena. We know of no other way to unite the principle of gauge symmetry with the 
description of massive vector bosons. While gauge symmetry is required in the theoretical 
description of electroweak physics in order to get meaningful predictions at the quantum 
level, the existence of the massive vector bosons W and Z is an experimental fact since 
their discovery at LEP. Although the Standard Model (SM) of particle physics is well 
tested even at the quantum level in several cases, there is no sign of the Higgs boson so 
far. Up to now only bounds on the Higgs boson mass(es) could be extracted from high 
energy collision experiments. These bounds are always obtained under the hypothesis that 
a certain model of the Higgs sector describes the data. However, as there is no Higgs signal 
yet, the choice of model for the Higgs sector is rather unconstrained. Therefore, extensions 
of the SM Higgs sector have to be considered seriously in phenomenology. 

Such extended Higgs sectors would also be considered as low energy effective theories 
of new physics models beyond the SM, which are proposed to describe physics at higher 
energies than the electroweak scale. Supersymmetry is one of the examples for such new 
physics scenarios, which requires at least two scalar doublet fields in the Higgs sector. 
Apart from supersymmetry, there are varieties of new physics models most of which have 
been proposed during the last decade, such as extra dimension models [T], the little Higgs 
models [2] and others. The low energy effective theories for some of these models predict 
extended Higgs sectors. Non-minimal Higgs sectors are also considered in models with 
strongly interacting dynamics for the symmetry breaking such as top-color models [3]. 
Furthermore, some models designed to explain tiny neutrino masses jHj, CP- violation [1], 
and electroweak baryogenesis [H] also require an extension of the minimal Higgs sector of 
the SM. 

From the phenomenological point of view, such extensions of the Higgs sector have to 
meet two major restrictions from experiments: a) the electroweak rho-parameter has to be 
one up to a few per mille and b) large flavor changing neutral currents (FCNC) have to 
be absent. By discarding models which do not meet these criteria or meet them only by 
fine-tuned choices of model parameters, one ends up with models which have p = 1 and 
no FCNC at tree-level. In particular, it is well known that all models with an arbitrary 
number of isospin SU (2)-doublets and -singlets can be of that type [Sj by imposing discrete 
symmetries to avoid FCNC. 

The Two Higgs Doublet Model (THDM) is the model with the minimal extension of 
the SM Higgs sector, which leads to new phenomena in the matter and gauge field sector 4 . 
In the THDM, tree level FCNC can be eliminated by imposing a discrete symmetry, under 
which there are two possibilities depending on its charge assignment [3 El; i- e -> (Type I) 
only one of the Higgs doublets gives mass to the fermions, and (Type II) one gives mass to 
the up-type quarks and the other to the down-type quarks and charged leptons. The Higgs 
sector of the minimal supersymmetric extension of the SM (MSSM) is a special case of 

4 A mere addition of scalar isospin-singlets would only modify the Higgs boson self-interactions. 



1 



the type II THDM. As the general THDM may be considered as a low energy description 
of other new physics models, it is an interesting question whether one can distinguish 
experimentally between the MSSM and such a model. 

The question about the new physics model underlying the Higgs sector can be addressed 
only after the existence of Higgs particles has been established experimentally. The Higgs 
search program at high energy colliders divides basically into three steps: a) discovery 
of the Higgs boson(s), b) measurement of fundamental properties like its mass, width, 
spin, parity etc., and c) determination of the underlying model of the Higgs sector; i.e., 
the measurement of the Higgs bosons' couplings to other particles and to themselves and 
the measurement of quantum effects. The current paper deals with an example which 
might contribute to last step. As is well known, the discovery of a charged Higgs boson 
would be an unambiguous sign of an extended Higgs sector. At hadron colliders, the main 
production mechanisms are top-quark decay in top-quark pair production if m H ± < mt+rrib 
jH] and single H~(H + ) production by bottom-gluon scattering and gluon fusion (gb — > tH~ , 
gg — > H~tb/ 'H~fv T and charge conjugated) if m H ± > m t + nif, [TUj . 

In this paper, we study H ± W T production at hadron colliders. This process is not 
a main production process for the charged Higgs boson, but it turns out to be strongly 
model dependent in contrast to the main production mechanisms. After the discovery of a 
charged Higgs boson, the observation of this process could potentially help to unravel the 
underlying model of the Higgs sector. The main purpose of this paper is to demonstrate 
this possibility by studying the predictions of the MSSM and THDM in comparison. The 
MSSM prediction for H ± W T production at hadron colliders has been studied by several 
authors [HI [T21 H3 HU H3 CHI IHI Especially, the detectability of H^W* production at 
the CERN Large Hadron Collider (LHC) in the framework of the MSSM has been studied 
in Ref. which concludes with rather poor prospects. We reconsider this process in the 
THDM, and study discriminative features with respect to the MSSM. The cross section 
is evaluated in the THDM in a wide parameter range. We exclude areas of parameter 
space by the requirement of vacuum stability [THl HH] and perturbative unitarity [20J |2TJ 
|2*2*] . 5 Experimental results, such as on the rho-parameter j^ZI, the muon anomalous 
magnetic moment [2H| and b —>■ sj [21] are also taken into account. We find that in the 
regions of parameter space which are not excluded by the theoretical requirements and 
the experimental constraints, a large enhancement of the hadronic cross section can be 
obtained as compared to the MSSM prediction. Therefore, in certain cases, the signal can 
be detectable at the LHC. 

The rest of the paper is organized as follows. In section 2, we briefly review the 
THDM and present the restrictions for the THDM parameters, derived from theoretical 
constraints and present experimental data. Section 3 deals with the essentials of the W ± H Zf 
production process at hadron colliders. In Section 4, we discuss the THDM prediction for 
the production cross section in comparison with the MSSM prediction, having regard to 
all the parameter restrictions discussed in section 2. Our conclusions and the appendix 
follow. 

5 Similar approaches are seen in Refs.|23| 
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2 THDM parameters 



The THDM contains two scalar weak isospin doublets $1, $2 with hypercharge Y(<&i) = +1. 
Conventionally, there are two types of THDMs which are characterized by the way the 
scalar doublets are coupled to fermions. In type I all fermions couple to one doublet, while 
in type II up- and down- type fermions couple to different doublets |H] • In either type FCNC 
are automatically absent at tree level j7j. Here, we consider a THDM of type II. The Higgs 
sector of the CP-conserving THDM can be described by the following potential [B]: 

V{$i, $ 2 ) = ^($1$! - vl) 2 + A 2 ($ 2 $ 2 - vlf + A 3 - vl) + ($ 2 $ 2 - vl)] 2 

+ A 4 [($ t 1 $i)($ 2 $ 2 ) - ($ t 1 $ 2 )($ 2 $ 1 )] + A 5 [Re($i$ 2 ) -v lV2 ] 2 + A 6 [lm($ t 1 $ 2 )] 2 , (1) 

with two real parameters with mass dimension one (f 1, f 2 ) and 6 dimensionless real param- 
eters (A 1; • • • , A 6 ). This potential has its minimum at $j = (0, Vi) T {i = 1,2). Re-expressing 
the fields such that the new degrees of freedom vanish at the minimum and diagonalizing 
the resulting bilinear scalar interaction terms one obtains a separation of the physical and 
unphysical spectrum of the model. The physical spectrum of the Higgs sector of the CP- 
conserving THDM consists of 3 neutral Higgs bosons, 2 CP-even (h°, H°) and one CP-odd 
(-4°), and two charged ones (H + , H~). The 8 free parameters of the Higgs sector can also 

be chosen to be the modulus of the vacuum expectation value v = \jv\-\-v\ = J 2~ 3 / 2 Gp 1 , 
the masses of the Higgs bosons, m^o, m#o, m^o and m#±, the mixing angles a and f3 used 
in the diagonalization of the Higgs boson propagator matrices, and A5. The discrete sym- 
metry of the Higgs potential (e.g. $1 — > $1, $ 2 — > — $ 2 ) is broken by a mass- dimension 
two term proportional to M 2 = v 2 X^. The physical meaning of M is the cut-off scale of 
the effective SM when M ^> v. In the following, we review the constraints which we take 
into account in our numerical study. 

1. theoretical constraints 

The coupling constants in Eq. (JTJ can be restricted by imposing theoretical requirements 
for the consistency of the model. Here, we use the conditions derived from the requirement 
of vacuum stability [THIHH and perturbative unitarity [20112111221 for the tree-level coupling 
constants. The condition of vacuum stability is given by 

Ai + A 3 > , A 2 + A 3 > , 

2 v / (Ai + A 3 )(A 2 + A 3 ) + 2A 3 + A 4 + min [0, A 5 - A 4 , A 6 - A 4 ] > . (2) 

The requirement of perturbative unitarity demands that the magnitudes of all tree-level 
S-wave amplitudes for the elastic scattering of longitudinally polarized gauge and Higgs 
bosons stay within the limit set by unitarity. In our analysis, we consider the 14 neu- 
tral channels [21]. The expressions for the eigenvalues of the scattering matrix, a» (i = 
1, . . . , 14), are summarized in Appendix [BJ 

2. rho-parameter constraint 

The electroweak rho-parameter, which is one at tree-level in the SM, is a scheme-dependent 
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quantity beyond the leading order in perturbation theory . For models beyond the SM the 
definition p = PoPsm = m w/( m2 z c w) * s used in Ref. |2Zj in the MS scheme, where psm 
absorbs all SM radiative corrections and po parameterizes the extra new physics contribu- 
tions, i.e. po = 1 if the new physics contributions vanish. A value for po has been obtained 
from a global fit to electroweak precision observables [27] . 

po = 0.9998l°;°o?o , (3) 

where the error bar given above is at the 2a level. The only difference between the SM and 
the THDM is the Higgs sector. Thus, we can decompose the contributions to p(= Pthdm) 
into all contributions containing at least one virtual Higgs boson, Pthdm, Higgs; and all others, 
which coincide with the SM contributions: Pthdm = PsM,no Higgs + Pthdm, Higgs- Likewise, 
the SM rho-parameter decomposes into psm = PsM.no Higgs + PsM,Hi gg s- Hence, the deviation 
5po = po — 1 at the one-loop level is just the difference 5po = Pthdm, Higgs — PsM,Hi gg s- In 
our study, we constrain the THDM parameters such that the one-loop prediction for Spo 
stays within the range indicated by Eq. Q; i.e., —0.0012 < 5p < 0.0023. The specific 
formulas for this difference can be found in Refs. [""] 130] . 

3. constraint from a M 

The latest results on the measurement of the anomalous magnetic moment of the muon 
"2H], a M , suggest that the difference between measurement and SM prediction is [2*7] : 

Aa M := a^ xp - - af 1 = (25.7 ± 8.54 ± 4.90) • HT 10 , (4) 

where the first error is the total experimental error and the second is the theoretical 
uncertainty of the SM prediction. In the THDM, radiative corrections to processes without 
any external Higgs bosons always split into contributions without any virtual Higgs boson 
and others which involve at least one virtual Higgs boson. Thus, the difference between 
the THDM and the SM prediction for a M is given by the difference in the Higgs sector 
contribution, Sa^ = a™ DM ' Hlggs — a SM,Higgs_ - g we jj known, the leading order virtual 
Higgs contributions are suppressed by two Higgs- muon Yukawa couplings [34J. At the 
two-loop level, a virtual Higgs boson can couple to the muon and one internal loop of 
heavy particles, e.g. Fermions with a much larger Yukawa coupling. This is the class 
of so-called "Barr-Zee" type Feynman graphs [S3]. Their contributions to a M exceed the 
leading order virtual Higgs contributions by several orders of magnitude. In our calculation 
of Sdfi, we include all relevant one-loop contributions [23] and all Barr-Zee type two-loop 
contributions with a closed fermion or charged Higgs boson loop [2""]. Assuming the validity 
of the THDM, we constrain its parameter space such that the value of 5a^ stays close to 
Aa M within the 2a error bars; i.e., -1.2 • 10" 10 < 5a M < 52.6 • 10~ 10 . 

4. constraint from b — > S7 

In our discussion, we take a charged Higgs boson with m H ± = 400 GeV. Therefore, the 
limits for new physics contributions to the decay b — *■ S7 are respected [3 [22] for all values 
of tan (3 which we discuss. 
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Figure 1: Typical Feynman graphs for the partonic processes gluon fusion and bb annihi- 
lation contributing to W ± H^ production at a hadron collider. For gluon fusion there are 
in total two triangle-type topologies for each quark flavor and six box-type topologies for 
each quark generation. 

3 W ± H T production at hadron colliders 

Partonic processes 

In the framework of the parton model, there are two distinct subprocesses which contribute 
to the production of a charged Higgs boson, H + , and a electroweak gauge boson, W: bb- 
annihilation, 

b(k, a, a) + b(k, (3, a) -> W~(p, A) + H + (p) , (5) 

and gluon fusion, 

g{k, a, a) + g(k, b, a) -> W~ (p, A) + H+(p) . (6) 

In- and outgoing momenta of the initial and final state particles are denoted by k, k and 
p, p respectively, the helicities of the initial state gluons or quarks by cr, a and of the 
final state W by A. a, b denote the gluon SU(3)-color indices and a, (3 the quark color 
indices. The square of the center of mass energy of the parton system is then given by 
s = (k + k) 2 = (p + p) 2 . 

The leading order ^-annihilation amplitude (0) consists of two types of tree-level Feyn- 
man graphs: (i) graphs with s-channel exchange of a neutral Higgs boson and (ii) graphs 
with a t-channel virtual top quark (see Fig. [TJ. The amplitude for the gluon-fusion process 
(JHJ) is given in leading order in perturbation theory by a set of one-loop Feynman graphs 
(see Fig|T]). Gluon fusion, though loop- induced, may contribute significantly to the cross 
section, because of the large number of gluon-gluon collisions with sufficient center-of-mass 
energy to exceed the production threshold at high energy hadron colliders. 

Associated W production via bb annihilation and the gluon fusion have been dis- 
cussed at first for a THDM with MSSM parameter values including the loop contributions 
from top and bottom quarks in the approximation = [IT]. Therefore, this study only 
covered MSSM scenarios with small tan/3 and scalar quark masses which are sufficiently 
heavy to decouple from the loop contributions. This work has been extended by including 



5 



squark-loop contributions and a non-zero 5-quark mass |T21 ll3|. thus allowing the investiga- 
tion of the process for arbitrary scalar quark masses and values of tan/3. The gluon-fusion 
channel in the MSSM has also been studied in Ref. ^1] and for the ^-annihilation channel 
the super symmetric electroweak corrections ^Hj and the QCD corrections ^B] at one-loop 
order are known. A phenomenological study of the signal-to-background ratio for the semi- 
leptonic signature W ± H Z¥ — > W^tb — > bbW + W~ — > bbjjl + missing energy /momentum for 
a MSSM-like THDM (neglecting the box-loop contributions) has been performed in Ref. 
|17j . The authors of Ref. 07] showed that the background rate, which mainly comes 
from it production, overwhelms the signal rate by two to three orders of magnitude. A 
generic feature of the MSSM gluon-fusion amplitude is a strong negative interference be- 
tween triangle- and box-type quark-loop Feynman graphs. This behavior, already noted in 
Ref. [TT], is due to the relations between Higgs masses in the MSSM, and leads to a much 
smaller cross section than one would obtain by squaring triangle- or box-type quark-loops 
alone. 

The W ± H T production cross section at hadron colliders in the framework of the gen- 
eral THDM has not been studied in the literature 6 . We discuss this cross section in the 
framework of a general CP-conserving type II THDM, which contains the MSSM Higgs 
sector as a special case and therfore allows for straightforward comparisons between the 
models. Dropping the superpartner contributions and allowing the Higgs boson masses 
(to^o, m H o, m^o, rriH±) and the mixing angle in the Higgs sector (a) to be free param- 
eters, one obtains the expression for the partonic cross section cr gg ^ w - H + in the THDM 
from the results of Ref. [T3J. The formulas for a b i^ w -H+ m the THDM can be taken over 
e.g. from the MSSM calculation in Ref. As all Higgs boson masses in the THDM 

are free parameters, it may occur that their values lie above the threshold + rriH±- 
This gives rise to resonant s-channel contributions. Therefore, we calculate the widths T$ 
($ = h°, H°, A ) of the neutral Higgs bosons in leading order approximation and intro- 
duce them in the s-channel propagators of the Feynman graphs of the bb- annihilation and 
gluon-fusion process . Specifically, we make the following replacements for the propagator 
terms of neutral Higgs bosons in the formulas of Refs. ^21 EI] : 



In the THDM, there are parameter combinations in which one Higgs boson can decay into 
two others which gives a large contribution to its width. We take this into account by 
including all self-interaction contributions to the total width of the neutral Higgs bosons. 
The partial widths of neutral Higgs bosons associated with decays into SM particles are 
well known and have been taken over from Ref. [B] and the ones associated with decays 
into Higgs particles are given by 



1 



1 




s — m\ + im$T<$> 




m\ + m\ + nij — 2m\mf — 2mfmj — 2mfm| 





6 A preliminary study of the THDM gluon-fusion process has been presented in Ref. |26|. 
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with the Higgs self-couplings g , $ 1 * < $ . listed in appendix 1X1 

We describe above how to obtain the expressions for the partonic cross sections from the 
formulas of Refs. [T2"l ll3|. However, our calculation of the partonic cross sections has been 
performed independently with the help of the computer programs FeynArts and FormCalc 

Hadronic cross section 

The hadronic inclusive cross section for W~ H + production in proton-proton collisions at 
a total hadronic center of mass energy can be written convolution E21, 

v P p-*w-H++x = y~] [ dr^f^ o nm _+ w - H +{TS,ot s {HR)) = Y] [ dVI^^ , (7) 

{n,m\ u {n,m\ v u 



with the parton luminosity 
dCPJL f 1 dx 1 



dr /_ x 1 + 5 n 



T T 
fn/p(x, flF)fm/p(-, ^F) + fm/p( x , ^F)fn/p(-, ^f) 



where f n / p (x, hf) denotes the density of partons of type n in the proton carrying a fraction 
x of the proton momentum at the scale /i^. In our case, there are two parton subprocesses 
contributing to inclusive W~H + hadroproduction; gluon fusion and ^-annihilation. The 
numerical evaluation has been carried out with the leading order MRST parton distribution 
functions [33] and with the renormalization and factorization scale fiR and fiF chosen equal 
to the threshold, m w + m#± . 



4 Numerical Results 

We take m Zl m w and G F as the input electroweak parameters, and use values mz = 
91.1876 GeV, m w = 80.424 GeV and G F = 1.16637 x 1(T 5 GeV" 2 [2Zj- For the strong 
coupling constant as([iR), we use the formula including the two-loop QCD corrections for 
rif = 5 with Aq CD = 174 MeV which can be found in [27j. The mass of the top and bottom 
quarks are fixed here as m t = 174.3 GeV and = 4.7 GeV. 

In our exemplary discussion we take one MSSM scenario as a reference, characterized 
by the following settings of MSSM parameters: 

Msusy = Mq = M = M D = M L = M E = 1000 GeV , 

[i = 300 GeV , (9) 
X t = A t - ii cot /3 = -1000 GeV , 
A u = Ad = A c = A s = Ah = , 

where Msusy is a common squark mass scale, /i is the Higgs-superfield mass parameter in 
the superpotential, X t the mixing parameter in the stop sector, and the A q are trilinear 
Higgs couplings to squarks 7 . We set m H ± = 400 GeV throughout the paper and calculate 

7 For more details on these parameters in our convention see Ref. |13| . 
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tan/? 






1.5 


3 


6 


10 


rriyfi [GeV] 


85.6 


105.5 


115.2 


117.6 


m H o [GeV] 


404.7 


396.8 


393.3 


392.4 


m A o [GeV] 


391.8 


391.8 


391.8 


391.8 


m H ± [GeV] 


400.0 


400.0 


400.0 


400.0 


a 


rad] 


-0.63465 


-0.36335 - 


0.19050 


-0.11558 


MSSM rrul 


37.78 


9.48 


3.08 


2.96 



Table 1: MSSM values for the Higgs masses and the mixing angle a. The prediction for 
the hadronic cross section in our MSSM reference scenario (see Eq. Q) are also shown. 



the mass of the MSSM CP-odd Higgs boson using the tree- level relation m 2 A0 = m 2 H± — m^. 
A THDM with all parameters set to the MSSM values, especially M = ttiao, is called 
MSSM-like. In Table we list the MSSM values for the Higgs masses and the mixing 
angle a. Note that the case tan/3 = 1.5 in the MSSM is actually ruled out by the latest 
LEP direct search results |2*H] . 

Fig. |21 shows a comparison between the THDM and MSSM predictions for the differ- 
ential hadronic cross section 

a nm ^ w - H +(VS, a s (fx R )) , (10) 

as a function of In the linear plot on the right hand side, one can see that the gg- 
and 66-channel contribution to the hadronic cross section is of comparable size, whereas 
the partonic cross section a gg _ rW - H + is orders of magnitude lower than a b i^ w - H +. This 
is caused by the enhancement factor due to the large number of gluon-gluon collisions at 
high energy hadron colliders. In Fig. |21 two THDM scenarios are displayed: one (scenario 
A) is completely MSSM-like, with the settings m H ± = 400 GeV and tan/3 = 6, and the the 
other (scenario B) coincides with the first one except for the choice m^o = 4m^ SSM , which 
leads to the peak at vl = m^o ~ 1600 GeV. We also show in Fig. El the MSSM prediction 
by thin solid and dashed lines further marked by circles and boxes respectively. Clearly, 
there is almost no deviation from the MSSM-like THDM except for squark threshold effects 
around \fl ~ 2000 GeV which can be seen in the left plot. 

In order to demonstrate the strong negative interference in the MSSM gluon-fusion 
process, we display artificial MSSM predictions for this process using only Feynman graphs 
either of box-type (thin dotted lines) or of triangle-type (thin dot-dashed lines) in the 
calculation. Those are more than one order of magnitude larger than the full result over 
a wide range of vl (see Fig. 121 left plot). From Fig. El we learn that we can get a 
much larger cross section near to the production-threshold region by near-decoupling of 
the triangle-type graphs. In this region the cross section in the THDM scenario B is close 
to the "only box graphs" MSSM result. However, the resonant peak far off the production 
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Figure 2: The differential hadronic cross section da nm /d\/§ in fb/GeV for W~H + produc- 
tion is plotted versus V§. Thick and thin lines show THDM and MSSM results respectively. 
Dashed and solid lines show the gg and bb contribution respectively. Two THDM scenarios 
are displayed: one with all parameters MSSM-like (A) and one with rriA = 4m^ SSM (B). 
The MSSM prediction is highlighted by circles (bb) and squares (gg). Also shown is the 
MSSM "prediction" for the gg cross section only using box- (dotted lines) or triangle-type 
(dot-dashed) Feynman graphs in the calculation. 



threshold also contributes significantly to the total hadronic cross section. Naively, one 
could have thought the peak far off would not contribute significantly because of the much 
lower parton luminosity. 

The hadronic cross section a pp ^ w ± H ^ as a function of m A o and m H o, for tan/3 = 1.5, 3, 
6, 10 is displayed for M 2 = m 2 A0 in Fig. El M 2 = m 2 A0 /2 in Fig. H and M 2 = in Fig. In 
all Figures the exclusion limits (thick lines) are superimposed on the cross section contours 
(thin lines). As a generic feature, the cross section rises strongly, if m^o or m#o becomes 
larger than mw+fnH ± ~ 480 GeV, with a maximum well above the contour where resonant 
A - and if °-exchange contributions are possible. This behavior is due to two enhancement 
effects: a) resonant propagator contributions for m A o,m H o > 480 GeV and b) a reduction 
of the strong negative interference if one propagator starts to decouple from the production 
threshold. The variation of the cross section with M is due to the M-dependence of the 
Higgs self-coupling constants (see Appendix EI), which enter the widths of A and H°. 

Imposing the constraints described in the previous chapter, we focus on the remaining 
areas in Figs. El to which are allowed by all constraints. Remarkably, there remains an 
allowed region in all cases. For M 2 = m 2 A0 (see Fig. EJ) the cross section in the allowed 
region varies roughly in the range 20 to 1200 fb for tan/? = 1.5, 10 to 50 fb for tan/3 = 3 
and around the MSSM value for tan/3 = 6 and 10. For M 2 = m 2 A0 /2 (see Fig. 0J the 
lower bound of the cross section is similar to the previous case, while 1000 fb are possible 
for tan/3 = 1.5, 3 and 6, and 20 fb for tan/3 = 10, which is still about 60 times the MSSM 
cross section. Depending on tan /3, the case M 2 = is strongly constrained to rather low 



9 



values of m#o. Because of the rho-parameter constraint, the allowed areas for tan/3 = 3, 6 
and 10 have a cross section below the corresponding MSSM scenario but of the same order. 
For tan/3 = 1.5 again large variations of the cross section are possible, roughly between 40 
and 1500 fb. 

5 Conclusions 

We have discussed the H ± W^ production at hadron colliders in the framework of a general 
type II THDM and compared the predictions with the MSSM. We find that the THDM 
prediction for the hadronic cross section can be completely different from the MSSM pre- 
diction. Specifically, we find regions in parameter space where the cross section exceeds 
1000 fb. These regions of large cross section are in agreement with theoretical constraints, 
vacuum stability and perturbative unitarity, and are not excluded by experimental con- 
straints from measurements of the rho-parameter and the muon magnetic moment. In 
Ref. [T7j, it is discussed that the size of the MSSM cross section would not be sufficient to 
detect this process at the LHC An enhanced cross section in a general THDM scenario, 
examples of which are shown in this paper, would give a nice possibility to see the H ± W Zf 
signal at the LHC. Certainly, the H ± W T cross section is not a discovery channel for the 
charged Higgs at hadron colliders. However, once a charged Higgs boson has been discov- 
ered, its observation will help to gain information on the underlying model of the Higgs 
sector. One should keep in mind that there are many other new physics models, apart 
from supersymmetry, which can be described by a THDM as a low energy effective theory. 
It would be valuable to study the detectability of H ± W T production in the framework of 
the general THDM by realistic simulation. We provide a FORTRAN code for general use. 
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Appendix 



A Higgs self-couplings 

We display here only those triple-Higgs couplings which are needed in the evaluation of the 
corresponding partial decay widths of the neutral Higgs bosons. The case of a MSSM-like 
Higgs sector is recovered if the MSSM values are chosen for m^o, m.#o, a and M is set to 
m^o. The shorthands := smip,^ := cosip are used. 



9h°H°H° 
9h°A°A° 
9h°H+H- 
9H°h°h° 

g H °A°A° 
9h°h+h- 



Gf 
Gf 
G p 
G p 
G p 



{^{s«-p - 3^ + 3£*!H>) + (4< + 2m 2 h0 )(^ - 4^) 

I V 2 Cp 0/33(3 



C/3 



(M 2 



-<a+f3 



I y/2 CpSp 
I a/2 CpSp 

\M 2 



+ Am AO s a -f3 - 2m 



h° ( 



c a+f3 



+ S 



a-/3) {, 



+ Am H± s a _f3 - 2m h0 ( h s a _ f3/ 



Sa+P + 3- 



3s 



s 



CpSp 

) - ( 4 «V + 2 %o)(- 



2 ■So—/? ■ 

'a , 

SpCp 



f M 2 S a+ p 2 2 , S a+ (3, 

Am AO c a _p + 2m H0 (c a _p , 



r A^f^ Sc-i-fl .2 n 2 / 

- Am H± c a -/3 + 2m H o{c a -/3 - 



CpSp 
S a+P ■ 

CpSp' 



B Higgs and gauge boson scattering amplitudes 

The tree-level S-wave amplitudes for the scattering of Higgs and longitudinal gauge bosons 
in the high energy regime can be calculated elegantly using the equivalence theorem |2lH l37|. 
We consider 14 distinct neutral channel scattering processes to be studied [2*T] . 

W£W£, W£H~, W£E + , E + E~,Z l Z l , Z l A, AA, Z L h, Z L E, Ah, AH, hh, hE, EE, 
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which lead to a 14-dimensional S-matrix, with the following eigenvalues, generically referred 
to as di(i = 1, . . . , 14), expressed in terms of the parameters Aj of the Higgs potential (JTJ). 



a± = T ^{3(A 1 + A 2 + 2A 3 ) ± (VH^i ~ A2) 2 + (4A 3 + A 4 + A 5 /2 + A 6 /2) s 

l07T L 

b± = + A 2 + 2A 3 ± V(Ai - A 2 ) 2 + (-2A 4 + A 5 + A 6 ) 2 /4} , 

c± = d ± = ^{A! + A 2 + 2A 3 ± ^(Ai - A 2 ) 2 + (A 5 - A 6 ) 2 /4} , 

L| 2 A3-A4-A 5 /2 + 5A 6 /2} , (12) 



61 = 16VT 



e 2 = y^{ 2A s + A 4 - A 5 /2 + A 6 /2} , 
/+ = Y^{ 2A 3 - ^4 + 5A 5 /2 - A 6 /2} 

/- = y^{ 2A3 + A * + As / 2 ~ V 2 } ' 

/i = / 2 = T ^{2A3 + A 5 /2 + A 6 /2}. 



In our study, we require the moduli of all these amplitudes to stay below 1/2 j|. 

To obtain the S-wave amplitudes ()12j) in terms of masses and mixing angles in the Higgs 
sector, we eliminate Aj(i = 1, 2, 3, 4, 6) according to the following relations. 



v 2 c 2 p 8v 2 V cj s 2 p ) 4 



A 2 — - 9 2 

A S 2a A5 m ff± \ _ m A° 

4f 2 s 2/ 3 4 ir v z 



with £ = m 2 H0 + m^o and A = m 2 H0 — m 2 ^. 
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Figure 3: Contour of the Hadronic cross section a pp ^w ± HT in fb for the case M 
(thin dotted lines). The limiting contours for 5po (thick solid lines) allow the cross-shaped 
area, and for 5a^ (thick dot-dashed lines) the whole displayed area for tan f3 = 1.5, 3 and the 
band-shaped area for tan/3 = 6, 10. The perturbative unitarity constraint, max |cij| < 1/2, 
(thick dashed lines) allows the pointed area which includes the lower left corner. The area 
allowed by the vacuum stability condition (0) is left of the thick short-dashed contour. The 
cross, labeled "MSSM", shows the point where all Higgs sector parameters coincide with 
the ones of our reference MSSM scenario. 
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Figure 4: Contour of the Hadronic cross section a pp ^w ± HT in fb for the case M 2 = m 2 A0 /2 
(thin dotted lines). The area of allowed parameter space, defined by the superposition 
of limiting contours for Spo (thick solid lines), Sa^ (thick dot-dashed lines), perturbative 
unitarity (thick dashed lines), and vacuum stability (thick short-dashed lines), is described 
in Fig. El 
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Figure 5: Contour of the Hadronic cross section a pp ^w ± HT hi fb for the case M 2 = (thin 
dotted lines). The limiting contours for the allowed values of Sp (thick solid lines) allows 
the cross-shaped area, for 5a M the whole displayed area for tan/3 = 1.5, 3 and the area 
below the thick dot-dashed line for tan/3 = 6, 10. The perturbative unitarity constraint, 
max|aj| < 1/2, (thick dashed lines) allows the almost rectangular area which includes the 
lower left corner. The vacuum stability condition (j2j) gives no constraint on the displayed 
area. 
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